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ABSTRACT 



We observed the low- mass X-ray binary and Z source GX 17+2 with the Rossi X- 
ray Timing Explorer during 6-8 Feb 1997, 1-4 Apr 1997, and 26-27 Jul 1997. The 
X-ray color-color diagram shows a clear Z track. Two simultaneous kHz quasi-periodic 
t-h I oscillations (QPOs) are present in each observation, whose frequencies are well correlated 

with the position of the source on the Z-track. At the left end of the horizontal branch 
(HB) only the higher frequency peak is observed, at 645±9 Hz, with an rms amplitude 
of 5.7±0.5%, and a FWHM of 183±35 Hz. When the source moves down the Z track 
to the upper normal branch the frequency of the kHz QPO increases to 1087±12 Hz, 
and the rms amplitude and FWHM decrease by a factor 2. Further down the Z the 
QPO becomes undetectable, with rms upper limits of typically 2.0%. Halfway down the 
HB a second QPO appears in the power spectra with a frequency of 480±23 Hz. The 
frequency of this QPO also increases when the source moves along the Z track, up to 
P H i 781±11 Hz halfway down the normal branch, while the rms amplitude and FWHM stay 

approximately constant at 2.5% and 70 Hz. The QPO frequency difference is constant 
at 293.5±7.5 Hz. Simultaneously with the kHz QPOs we detect HB QPOs (HBOs). 
The simultaneous presence of HBOs and kHz QPOs excludes the magnetospheric beat- 
frequency model as the explanation for at least one of these two phenomena. 
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1. Introduction 

GX 17+2 is a bright low-mass X-ray binary and 
Z source (Hasinger & van der Klis 1989). Along the 
Z track traced out in the X-ray color-color diagram 
(CD) the mass accretion rate (M) is thought to in- 
crease from the horizontal branch (HB), via the nor- 
mal branch (NB), to the flaring branch (FB). On the 
HB and upper NB 18-60 Hz quasi-periodic oscilla- 
tions (QPOs) are present (called HBOs), and on the 
lower NB and FB 5-20 Hz QPOs (N/FBOs). Atoll 
sources have no Z track in the CD (Hasinger & van 
der Klis 1989), and do not exhibit HBO or N/FBOg. 
However, many show kHz QPOs between 300 and 
1200 Hz (see van der Klis 1997 for a recent review). 
Usually two simultaneous kHz QPOs (hereafter: twin 
peaks) are seen, whose frequencies increase with M , 
while their difference (the twin peak separation) re- 
mains constant. In X-ray bursts oscillations are seen 
with a frequency once or twice the twin peak sepa- 
ration. This strongly suggest a model in which the 
burst oscillation frequency is the spin frequency (or 
twice that), the higher kHz QPO frequency some in- 
ner disk frequency, and the lower kHz QPO frequency 
their beat. 

So far, two Z sources have shown kHz QPOs: Sco 
X-l with frequencies in the range 550-1100 Hz (van 
der Klis et al. 1996a; 1997b), and GX 5-1 (325- 
896 Hz; van der Klis et al. 1996b). Both sources 
showed twin peaks. In GX 5—1 the QPO frequencies 
increased from the left end of the HB to the upper NB 
(van der Klis et al. 1996b), in Sco X-l from the upper 
NB to the FB (van der Klis et al. 1996a; 1997b). In 
view of the similarities between the kHz QPOs in Z 
and atoll sources, both are most likely caused by the 
same mechanism. The only difference so far is that in 
atoll sources the peak separation is constant, whereas 
in Sco X-l it decreases when the QPO frequencies 
increase. In this Letter, we report the discovery of 
twin kHz QPOs in the Z source GX 17+2 (see also 
van der Klis et al. 1997a). 

2. Observations and Analysis 

We observed GX 17+2 with the PCA onboard the 
Rossi X-ray Timing Explorer on 6-8 Feb, 1-4 Apr, 
and 26-27 Jul 1997. A total of 120.5 ksec of good data 



Although in several of them broad QPO-like peaks are found 
with frequency around 20 Hz (e.g. Hasinger & van der Klis 
1989; Wijnands & van der Klis 1997) 



were obtained. During Feb 8 (1.1 ksec) and a small 
part of Feb 7 (7.8 ksec) and Jul 27 (4.8 ksec) only 4 
of the 5 detectors were on. The source then covered 
parts of the Z also well-covered with 5 detectors. As 
adding the 4-detector data did not improve the results 
significantly, we only used the 5-detector data to have 
the most accurate CD possible. Data were collected 
with 16 s time resolution in 129 photon energy bands 
(energy range 2-60 keV), and simultaneously with 122 
/zs time resolution in 4 bands: 2-5 keV, 5-6.4 keV, 
6.4-8.6 keV, and 8.6-60 keV. 

We made power density spectra from the 122 [is 
data, using 16 s data segments. To measure the 
kHz QPOs we fitted the 96-4096 Hz power spectra 
with a function described by a constant, one or two 
Lorentzian peaks, and a broad sinusoid to represent 
the dead-time modified Poisson noise (Zhang et al. 
1995; Zhang 1995). The Very Large Event window 
(van der Klis et al. 1997b) was set to 55 /is, so that 
its effect on the Poisson noise was small enough that 
it could be absorbed by the broad sinusoid. To mea- 
sure the HBOs we fitted the 8-256 Hz power spectra 
with a constant, one or two Lorentzian peaks, and 
a power law representing the continuum. Differen- 
tial dead time (van der Klis 1989) was negligible. We 
determined errors using A\ 2 = 1.0 and upper limits 
using Ax 2 = 2.71, corresponding to a 95% confidence 
level. Upper limits on the kHz QPOs were determined 
using a fixed FWHM (100 Hz). In the presence of 
one kHz peak the upper limit on the other was de- 
termined by fixing its frequency at the frequency of 
the detected QPO plus or minus the mean peak sep- 
aration (depending on whether the lower or higher 
frequency QPO was detected, respectively). Upper 
limits on the HBO second harmonic were determined 
using a fixed FWHM (10 Hz) and a frequency twice 
that of the first harmonic. 

In the CD and the hardness-intensity diagram 
(HID) we used for the soft color the log of the 3.5-6.4 
keV/2.0-3.5 keV count rate ratio, for the hard color 
the log of the 9.7-16.0 keV/6.4-9.7 keV ratio, and as 
intensity the log of the count rates in the 2.0-16.0 keV 
band. 

We used the S z parameterization (Wijnands et al. 
1997a, and references therein) for measuring the posi- 
tion along the Z. The HB/NB vertex is at S z =1.0, the 
NB/FB vertex at S z =2.0. By using logarithmic val- 
ues for the colors, S z does not depend on the values of 
the colors but only on their variations (Wijnands et 
al. 1997a). We selected the power spectra according 
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to S z , and determined the average S z for each aver- 
age power spectrum. The S z error bars represent the 
standard deviation of the S z distributions. 

3. Results 

In April the source traced out a full Z in the CD 
and HID (Fig. [[]). In February the HB/NB vertex 
and the upper NB were covered, and in July the HB. 
The Feb data fall right on top of the Apr data, indi- 
cating that the Z did not move appreciably between 
February and April. The Jul data have higher count 
rate and higher soft color. Such shifts have been ob- 
served before in GX 17+2 by Kuulkers et al. (1997), 
and also in other Z sources (Cyg X-2: Hasinger et 
al. 1990; GX 5—1: Kuulkers et al. 1994). In all ob- 
servations we detected kHz QPOs in the 8.6-60 keV 
band. In the 6.4-8.6 keV band the QPOs were not, 
or marginally, detected. Combining both bands the 
significance of the QPOs increased, although the rms 
amplitude decreased a little. Therefore, we used the 
6.4-60 keV range for the analysis of the QPOs. 

We detected the kHz QPO only for S z < 1.5. Typ- 
ical power spectra are shown in Fig. ^. Note the si- 
multaneous presence of the HBO and its harmonic 
(upper frame), and the kHz QPOs. The proper- 
ties of the QPOs as a function of S z are shown 
in Fig. |^. The higher frequency QPOs increase 
in frequency from 645±9 Hz at S z =0.07±0.05 to 
1087±12 Hz at S Z =1.44±0.Q5 (Fig.|a), while their 
rms amplitude and FWHM decrease from 5.7±0.5% 
to 1.5±0.4% (Fig|c), and from 183±35 Hz to 21±20 
Hz. Above S z =1.5 the higher frequency QPO is un- 
detectable (upper limits typically 2.5% rms). For 
S z <0.5 the lower frequency QPO is not detected, with 
upper limits near 2.5% rms. From S z =0.56±0.04 to 
S z =1.34±0.05 the lower frequency QPO is detected at 
a frequency increasing from 480±23 Hz to 781±11 Hz 
(Fig.|a). The rms amplitude (Fig.|e) and FWHM 
of the lower frequency QPO are about constant be- 
tween 2-3% and 50-150 Hz. When both QPOs are 
detected the frequency difference is consistent with 
being constant at 293.5±7.5 Hz. 

The QPOs were strongest at high energy. We com- 
bined all data with S z <0.5 (S z =0.23±0.14), and with 
S z between 1.0 and 1.3 (S z =1.15±0.10) to detect the 
QPOs at lower energy. When S z >1.0 both QPOs are 
seen. Both get stronger with energy and both have 
about the same amplitude in each band (see Tab. 0). 
When S z <0.5 only the higher frequency QPO is seen. 



Apart from an overall increase in amplitude its energy 
dependence is approximately the same as for S z >1.0 
(see Tab. [t]). 

The kHz QPOs and HBO and its second harmonic 
occur simultaneously (see Fig. ^|). The HBO fre- 
quency increases from the left end of the HB onto the 
upper NB (Fig. |d). At S z =1.34±0.05 it reaches a 
maximum of 61±0.2 Hz, further down the NB it de- 
creases again. On the HB the rms amplitude of the 
fundamental decreases with S z ; on the NB it is about 
constant at 1.5-2% (Fig.|d). Above S z =1.6 it is not 
detected with upper limits of 2.0% rms. The rms am- 
plitude of the second harmonic decreases on the HB 
with S z (Fig. ^|f). A detailed analysis of the HBO and 
N/FBO will be reported by Homan et al. (1997). 

4. Discussion 

We detected two simultaneous kHz QPOs in GX 
17+2. After Sco X-l (van der Klis et al 1996a; 1997b) 
and GX 5—1 (van der Klis et al. 1996b) this is the 
third Z source displaying kHz QPOs^. The QPO fre- 
quency increases down the HB and onto the NB. For 
atoll sources an increase in frequency is thought to 
be due to an increase in M. Due to the very simi- 
lar properties of the kHz QPOs in atoll sources and 
in GX 17+2 our observations confirm that M indeed 
increases when the source moves down the Z, even 
when (on the NB) the count rate decreases, as previ- 
ously concluded by Hasinger & van der Klis (1989). 

Although the behavior of the X-ray spectra changes 
significantly at S z =l when the source enters the NB, 
and the X-ray flux begins to decrease with S z there, 
the frequencies of all QPOs continue to increase with 
S z in the same manner as on the HB. The only ev- 
idence that the rapid variability is affected as the 
source enters the NB is by the HBO rms amplitude. 
Decreasing with S z on the HB, the amplitude is about 
constant on the NB. 

Between S z =1.4 and 1.5 the HBO frequency de- 
creases. A similar decrease in HBO frequency on the 
NB for GX 17+2 was seen in the "variable-frequency" 
QPO on the NB reported by Wijnands et al. (1996). 
Although supporting evidence for a decreasing HBO 
frequency on the NB was found in Ginga data on 
Cygnus X-2 (Wijnands et al. 1997a), it was unclear 
if the "variable-frequency" QPO in GX17+2 was in- 

9 Recently, we also found kHz QPOs in Cygnus X-2 (Wijnands 
et al. 1997b). 
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deed the HBO (see Kuulkers et al. 1997). From Fig. 

it is now evident that indeed it was, and also that 
it was the fundamental and not, as proposed by Ku- 
ulkers et al. (1997), the second harmonic. The data 
reported here and by Wijnands et al. (1996) show 
that the behavior of the HBO did not change signifi- 
cantly over one year. In the data presented here the 
HBO does not decrease to the same frequency as in 
1996; this is due to a lack of data on the NB. 

Beat-frequency models (BFMs) are the best mod- 
els so far to explain the kHz QPOs in atoll sources. 
In GX 17+2 kHz QPOs and HBOs are observed si- 
multaneously, so if the magnetospheric BFM explains 
the HBO (Alpar & Shaham 1985; Lamb et al. 1985) 
it can not explain the kHz QPOs (the same is true 
for Sco X-l; van der Klis et al. 1997b). The sonic- 
point BFM for the kHz QPOs in atoll sources (Miller, 
Lamb, & Psaltis 1997) assumes that the highest fre- 
quency is due to orbital motion around the neutron 
star at the sonic radius, and the lower frequency to 
its beat with the neutron star spin frequency. This 
model predicts a constant twin peak separation, as 
observed for atoll sources. According to this model 
the spin period of the neutron star in GX 17+2 would 
be 3.41+0.03 msec. 

In Sco X-l the twin peak separation decreases 
when the frequencies of the kHz QPOs increase (van 
der Klis et al. 1997b). Although the peak separation 
in GX 17+2 is consistent with being constant, its er- 
ror bars are larger than in Sco X-l and a decrease 
similar to that in Sco X-l can not be excluded (see 
Fig. [|) . Note that the decrease of the peak separation 
in Sco X-l was observed when this source was moving 
down the NB onto the beginning of the FB. We could 
not observe the kHz QPOs this far down the Z in GX 
17+2. 

Another explanation that has been proposed for 
kHz QPOs is the Photon Bubble Oscillation (PBO) 
model of Klein et al. (1996). These oscillations 
have been modeled to occur below 0.4 L Edd and only 
weakly (or not at all) above 1.1 Lgdd (Klein et al. 
1996; R. Klein, private communication 1997). The 
amplitudes of the PBOs decrease and their frequen- 
cies increase with M, just as we observe. However, in 
calculations so far PBOs near 1200 Hz occur at 0.4 
LEdd, a factor two lower than the inferred luminos- 
ity when we observed 1200 Hz QPOs. Calculations 
for luminosities between 0.4 and 1.1 LEdd covering 
a wider range of parameters are needed to compare 
the properties of the PBOs more accurately with the 



observed properties of the kHz QPOs in Z sources. 

GX 17+2 and Sco X-l have similar timing proper- 
ties. The highest observed frequency of the kHz QPOs 
in both sources is approximately 1100 Hz, and the 
peak separations are more or less equal (see Fig.|J). 
However, several aspects of the timing behavior differ 
significantly between the two sources. In Sco X-l the 
highest frequency for the kHz QPOs is observed at the 
beginning of the FB (van der Klis et al. 1996a), but in 
GX 17+2 this frequency is already reached on the up- 
per NB. This difference is remarkable, because both 
sources are thought to reach LEdd near the NB/FB 
vertex (i.e. S z =2). So, the same frequency is reached 
in GX 17+2 further from the Eddington critical value 
than in Sco X-l. 

Also, in Sco X-l the frequency of the HBO does not 
decrease as M increases (van der Klis et al. 1997b). 
Wijnands et al. (1996) suggested that an increase in 
the mass flux through the radial inflow as GX 17+2 
moves down the NB could produce the decrease of 
mass flux through the inner disk needed in the mag- 
netospheric BFM to explain the observed decrease of 
the HBO frequency with increasing M. It remains 
to be seen if indications can be found for a differ- 
ent radial inflow behavior between the two sources, 
as would be required to explain the difference in the 
HBO frequency changes with M within this frame- 
work. 
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Fig. 1. — Color-color diagram (a) and hardness- 
intensity diagram (b) of GX 17+2. The soft color 
is the logarithm of the ratio of the count rates be- 
tween 3.5-6.4 keV and 2.0-3.5 keV; the hard color is 
the logarithm of the ratio of the count rates between 
9.7-16.0 keV and 6.4-9.7 keV; the intensity is the log- 
arithm of the count rate in the photon energy range 
2.0-16.0 keV. The background, which is almost negli- 
gible (50 ct/s in the energy range 2.0-16.0 keV), was 
subtracted, but the intensity was not dead-time cor- 
rected. All points are 16 s averages. The thick line in 
the CD is the spline which was used to determine S z . 
Typical error bars are shown at one of the points. 



6 



CO 
Cxi 



CN - 



in 

Csi 



CM 



o 



CM 



m 
cn 




= 0.140±0.041 




I II 4- | wH I I I 1 1 1 1 1 1 I I I III 1 I I I 1 1 




S z =1 .059±0.044 




10 100 1000 

Frequency (Hz) 

Fig. 2. — Typical Leahy normalized power spectra 
in the energy range 6.4-60 keV. In the top panel, 
for S z =0.140±0.041, clearly the HBO, its second har- 
monic, and the higher frequency QPO at 675 Hz are 
visible. In the bottom panel, for S z = 1.059±0.044, the 
HBO fundamental, the lower frequency QPO (680 Hz) 
and the higher frequency QPO (986 Hz) are visible. 
The upward trend of the power at kHz frequencies, 
best seen in the bottom panel, is due to instrumental 
dead time. 
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Fig. 3. — (a) Frequency of the kHz QPOs, (b) fre- 
quency of the first and second harmonic of the HBO, 
(c) the rms amplitude of the higher frequency kHz 
QPO, (d) the rms amplitude of the first harmonic 
of the HBO, (e) the rms amplitude of the lower fre- 
quency kHz QPO, and (/) the rms amplitude of the 
second harmonic of the HBO as a function of S z . 
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Fig. 4. — The kHz peak separation versus the fre- 
quency of the higher frequency QPO. The bold points 
with the large error bars are for GX 17+2, the fine 
points with the small error bars are for Sco X-l (taken 
from van der Klis et al. 1997b, their Figure 3a). 
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Table 1 

kHz QPO rms amplitude versus photon energy 





S z =0.5±0.14 


S z =1.15±0.10 




Upper peak 


Upper peak 


Lower peak 


Energy 


rms 


rms 


rms 


(keV) 


(%) 


(%) 


(%) 


2-5.0 


< 1.9 


< 1.1 


< 1.2 


5.0-6.4 


3.4±0.5 


2.0±0.5 


< 1.9 


6.4-8.6 


3.9±0.5 


2.2±0.4 


1.7±0.5 


8.6-60 


7.4±0.2 


3.6±0.3 


4.1±0.2 


5.0-60" 


5.2±0.2 


2.1±0.2 


2.3±0.2 



a All errors correspond to A% 2 = 1. The upper limits 
correspond to the 95% confidence level. 

b Addcd for comparison with Cygnus X-2 (Wijnands 
et al. 1997b). 
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